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ABSTRACT 

Context. There is a continued debate as to the form of the outer disc of the Milky Way galaxy, 
which has important implications for its formation. Stars are known to exist at a galacto-centric 
distance of at least 20 kpc. However, there is much debate as to whether these stars can be 
explained as being part of the disc or whether another extra galactic structure, the so called 
Monoceros ring/stream, is required 
Aims. To examine the outer disc of the Galaxy toward the anti-centre to determine whether the 
star counts can be explained by the thin and thick discs alone. 

Methods. Using Sloan star counts and extracting the late F and early G dwarfs it is possible to 
directly determine the density of stars out to a galacto-centric distance of about 25 kpc. These are 
then compared with a simple flared disc model. 

Results. A flared disc model is shown to reproduce the counts along the line of sights examined, 
if the thick disc does not have a sharp cut off. The flare starts at a Galacto-centric radius of 16 kpc 
and has a scale length of 4.5 ±1.5 kpc 

Conclusions. Whilst the interpretation of the counts in terms of a ring/stream cannot be defi- 
nitely discounted, it does not appear to be necessary, at least along the lines of sight examined 
towards the anti centre. 
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1. Introduction 

The position of the Sun within the Galactic plane has made understanding the true form of the 
Galaxy difficult. In the last ten years the situation has improved, however, with the advent of deep 
large areas surveys. Surveys in the near IR such as 2MASS, UKIDSS or in the visible such as the 
Sloan Digital Sky Survey (SDSS, Newberg et al. 2002) have allowed a more precise modeling of 
the large scale structures of our Galaxy. 

It is generally accepted that there is both a thin and thick stellar disc. The stellar density of the 
thin disc is exponential in height above the Galactic plane and radial distance. In the vicinity of the 
Sun it has a scale length of about 2.4 kpc, however the scale height depends on the source type. The 
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older sources have a scale height of around 185 pc, but for the younger sources it is considerably 
less than this. The density of the thick disc (Gilmore & Reid 1983) is also exponential but with a 
larger scale length and scale height (Bilir et al. 2008), although it has less than a tenth of the density 
of the thin disc on the plane in the vicinity of the Sun. 

Of particular importance to our understanding of the Galaxy and its formation, is to determine 
what happens towards the outer edge of the disc. Robin et al. (1992), using visible images from the 
CFHT, suggested that Galaxy has a cut of at galacto-centric radius (R) of 13.5 kpc when looking at 
/ = 179° b = -2.5°, as the stellar density drops rapidly beyond that point. Lopez-Corredoira et al. 
(2002) using 2MASS data suggested that the drop in stellar density near the plane was not caused 
by a cut off, rather a flare in the disc and concluded that the disc did not have a well defined cut off 
to at least 15 kpc. 

The most recent deep surveys such as the SDSS clearly show that there are stars out to at least 
R-2Q kpc. Their distribution, however, is not consistent with a simple exponential disc particularly 
when looking 15 to 30 degrees from the plane. Some authors (Momany et al. 2006) suggest these 
can be explained by a flare of the disc. Many other authors, however, have preferred to attribute 
these sources to rings or streams beyond the edge of the disc. In particular, the sources being 
studied here have been attributed to the so called Monoceros ring/stream (Newberg et al. 2002; 
Rocha-Pinto et al. 2003; Conn et al. 2005, 2007, 2008). This has been associated to the remnants 
of a dwarf galaxy which was cannibalized by the Milky Way, and whose progenitor was associated 
to the Canis Major over-density (Martin et al. 2004; Martmez-Delgado et al. 2005; Bellazzini et 
al. 2006; Butler et al. 2007; de Jong et al. 2007; Conn et al. 2007). This ring runs approximately 
parallel to the galactic plane, in the latitude range 10 < \b\ < 35 deg. and over most of the second 
and third quadrants. 

That stellar streams exist is well known, and that produced by the Sagittarius dwarf galaxy 
(Ibata et al. 1994) is probably the most well studied stream associated with the Milky Way. Its 
form and high angle to the Galactic Plane makes confusion with a Galactic component highly 
unlikely. However, the proposed Monoceros stream is a very different feature to the Sagittarius 
stream running, as it does, parallel to the Galactic plane. Furthermore, whilst the Sagittarius stream 
is clearly a dwarf galaxy, the progenitor of Monoceros in the over-density in Canis Major has also 
been questioned, as other authors have explained the excess star as an effect of the warped+flared 
disc of the Milky Way (Momany et al. 2004, 2006; Lopez-Corredoira 2006; Lopez-Corredoira et 
al. 2007). 

In this paper we examine if, using a simple flared disc model as suggested by Momany et al. 
(2006), we can reproduce the form of the deep star counts seen 15 to 30 degrees off the plane in 
Galactic anti-centre region, without requiring the presence of extra-galactic streams. 

2. The Method 

For this work we wish to count sources at a R < 25 kpc towards the Galactic anti-centre and 
between 1 1 and 3 1 degrees off the plane. 

By far the simplest method of determining the stellar density along a line of sight in the disc 
is by isolating a group of stars with the same colour and absolute magnitude within a colour mag- 
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nitude diagram. This allows the luminosity function to be replaced by a constant in the stellar 
statistics equation: 

. dN(m) In 10 3 

Mm) = — = ——ujp[r(m)]r(my, (1) 

dm 5 

r(m) = 10['"- M+5 ]/ 5 , 

where w is the area of the solid angle in radians and r is the distance in parsecs; the differential star 
counts for each line of sight, A(m), can be immediately converted into density p(r). 

In the near IR, the red clump giants have been successfully used as standard candles, partic- 
ularly when looking in the inner Galaxy (Lopez-Corredoira et al. 2002). However, the red clump 
stars in the outer disc at the distance of interest would appear at nik-15 and so the local dwarfs with 
the same colour would completely dominate the counts at this magnitude (Lopez-Corredoira et al. 
2002). 

As the areas of interest all lie well off the plane, then the extinction will be low (typically 
under 1 magnitude). Therefore, it is possible to use visible counts. An examination of the HR 
diagram shows that when the extinction is low, the late F and early G dwarfs can be isolated using 
colour with only minimal contamination from other sources with the same colour but different 
absolute magnitudes (sub dwarfs,galacto-centric radius giants etc). For this work we have selected 
the source between F8 V and G5 V. This gives a range of g - r of 0.36 to 0.49 and a range of absolute 
magnitudes M g -4.2 to 5.4 which makes the sources approximately m g ~ 20 at the distance of 
interest. By having a range of absolute magnitudes, then when the counts are converted to density 
vs distance there is some smoothing and this is not included in a model which assumes a single 
absolute magnitude. Although the selection used here is not as well defined as the red clump stars in 
the near IR, it still has a sufficiently small range in absolute magnitude that the smoothing has little 
effect and the above approximation remains valid. There will, however, be sufficient stars detected 
in the outer Galaxy to give meaningful statistics which would not be the case if a smaller range 
in absolute magnitudes were used. Earlier sources have not been included as these sources would 
belong to a younger population with a far smaller scale height, also the absolute magnitude changes 
far more rapidly with colour. Later sources would have significant giant contamination and again 
the absolute magnitude changes more rapidly with colour. 

The areas selected for this work lie between 10 and 30 degrees from the plane and are within 
40 degrees of the anti-centre. We have not extended the range of longitudes because of the Galactic 
warp. The warp means that, in effect the position of the Galactic plane varies with position in the 
Galaxy, and along the line of sight. Although the effect of warp is strongest towards the outer edge 
of the disc, where the plane can be over a kpc away from the expected position with no warp, it can 
be clearly seen in the counts within a few kpc of the Sun (Hammersley et al 1995, Lopez-Corredoira 
et al. 2002, Momany et al. 2006). The models of the warp are normally simple, based on tilted rings, 
and in general these do provide a good representation of the star counts. However, in the very outer 
disc the models are harder to test, and a small error in the model can lead to a significant error to 
the predicted position of the plane, and hence the predicted densities. Therefore, we have limited 
the regions used for this work to those where the effect of the warp should be small and so can be 
ignored. It should be noted, however, that whilst it will be possible to ascertain whether a flare can 
reproduce the star counts seen in the outer disc, it will not be possible to accurately determine all 
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Table 1. Lines of sight used. 
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Fig. 1. The CMD for the region / = 183°, b = 21 c 



of the parameters. When large scale deep star counts in the outer disc become available it will be 
possible to model the warp more accurately. Until then, however, great care must be taken when 
comparing models with counts in the outer disc. 

A second reason for limiting the longitude is that the distance to the feature increases rapidly 
with distance from the anti-centre, making the sources fainter and appear closer to the plane so that 
there is less contrast with the other disc sources. 



3. The Data 

The source densities in these regions will be very low requiring square degrees of sky to be covered 
to provide sufficient counts to give reasonable statistics. Therefore, the data has been taken from 
the SDSS release DR7 (Abazajian et al. 2009) with regions ~ 1 square degree. As the regions are 
all well off the plane it has been assumed that the extinction is local and so all of the magnitudes 
have been corrected for extinction using the Galactic extinction model of Shlegel et al (1998). 
Furthermore, the extinction is relatively small and a small residual error will not significantly affect 
the results, although it would make the sources appear at the wrong distance. Table 1 shows the 
position and assumed extinctions. The positions / = 183°, b = 21° and / = 220° where used by 
Newberg et al. (2002); and the position at / = 150° was used by Conn et al. (2005), to support their 
"Monoceros stream" hypothesis. 
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Fig. 2. The differential star counts with the error bars in five regions towards the anti centre. As 
well as the counts, the predictions of the Besangon model and a simple flared or non-flared disc 
model are also shown. 

4. Results 

Figure 1 shows the g vs. g - r colour magnitude diagram for the region I = 183°, b — 21° with the 
sources of interest marked. It can be seen they are in a reasonably uncluttered part of the diagram 
allowing them to be isolated. Furthermore, there are no 'holes' in the area of interest indicting that 
there is no areas of high extinction. Figure 2 shows the differential counts for the selected colour 
along with the Besanfon model (Robin et al. 2003) for these stars. This model is widely used and 
generally accurately predicts the counts over wide areas of the sky. Here we have only run the 
model for the stars between F9V and G4V. The model was run with no extinction and predicted 
counts in V which were then converted to g using the g — V colour for these stars. The Poissonian 
error bars from the counts in each magnitude bin is also shown. 

The counts have two clear peaks, one at about m g = 17 and a second at m g =20.7 with a 
minimum at m g = 19.5 (see Fig. 2). This is seen most clearly at b = 21°, and b = 31°, however 
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Fig. 3. The stellar density with R along five lines of sight towards the anti centre. As well as the 
counts the predictions of the Besancon model and a simple flared or non-flared disc model are also 
shown. 



even in the region / = 180°, b — 1 1° there is some evidence for excess counts. The Besangon model 
does give a good fit to the data for m g < 18, however, at m g = 20 there is a sharp rise in the counts 
which is not predicted by the model. The peak corresponds with a distance of about 15 kpc from 
the Sun or a R = 23 kpc. This second peak cannot be caused by the halo, which is included in the 
model Besangon as the peak for these counts is predicted to start at m g « 23 and reach a peak at 
about m g = 25. Therefore, there are two possible alternatives; this is part of the disc or an extra 
galactic component, such as a stream beyond the disc of the Galaxy. 



5. Analysis 

These star counts cannot be produced by a simple exponential disc, even when including both the 
thin and thick disc, as such a model cannot produce this second bump (see figure 2). Even if there 
were no cut off this would only slightly affect the counts. However, a flared disc with no cut off 
may produce the second peak along specific lines of sight, as it would on average distribute stars 
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away from the plane. In order to explore this possibility, Eq. (1) has been used directly determine 
the stellar density with distance assuming that the luminosity function can be replaced by a delta 
function with M g = 4.80. Figure 3 shows the density with galacto-centric radius for the lines of 
sight used. It is notable that the stellar density is continuously decreasing, although it does change 
slope, without significant gaps or bumps suggesting that it may be explained only one continuous 
Galactic structure and not two separate structures. Hence the second bump in Fig. 2 is an effect of 
the shallower slope in the stellar density being multiplied by the factor r 3 [see Eq. (1)]. 

In order to determine the effect of a flare, a simple model was built containing a thin disc(scale 
height h z jhm,o = 186 pc and scale length h Rj Mn = 2400 pc) and a thick disc (scale height h a hick,o = 
63 1 pc and scale length h Rj Mck = 3500 pc) were assumed with the density of the thick disc being 
9% of that of the thin dick on the plane at the solar radius (as observed at high latitude with SDSS; 
Bilir et al. 2008). The halo counts of the model by Bilir et al. (2008) were also included, although 
they give a negligible contribution at these magnitudes. It was assumed that the flare increase the 
scale height of both the thin and thick discs beyond a certain galacto-centric radius, and that the size 
of the flare increased exponentially. It should be noted, however, that the thin disc only provides 
a significant contribution to the counts where the line of sight is within a few hundred pc of the 
plane. The flare begins at about R = 16kpc by which distance the lines of sight are over a kpc from 
the plane and hence the contribution of the thin disc to the counts where the flare is important is 
negligible. Therefore, this paper is only analysing the effect of the flare on the thick disc. Lines of 
sight far closer to the plane would be required in order to analyse the effect of the flare on the thin 
disc, and only then can any relationship between the flare in the thin and thick discs determined. 

It was also assumed that the total number of stars at a galacto-centric radius followed a simple 
exponential model. If this is not done then the large scale height caused by the flare would lead to 
there being more stars in total at that galacto-centric radius. 

The variation of stellar density is given by 

P = Pthin + P thick + Phalo, (2) 
A ( h z ,Mn,Q \ ( R-Rp\ ( Izl 

Pthin = A — exp — - exp - 



hz,thin(R)j \ h Rj hin ) \ h ZJ hin(R) 

Pthick = 0.09A 



exp 

Phalo = 1-4 x 10~ 3 A 



/ K Mck& \ I R-Rq \ x I jzj \ 

\h z ,thick(R) I \ flRjhick I \ h z ,thick(R)j 

Ha" 1 )] 



10.093 



(R sp /R Q yi* 

I h 7 jhi n lthick.Oi 

h z jhin/thick(R) = 



h z ,thin/thick,&, R < 16 kpc 

h z ,M n /thick,0 exp ( fi '^ kpc ) , R > 16 kpc 



R sp = Vfl 2 + 2.52z 2 

This is, therefore, the standard thin and thick exponential disc model, but with a flare starting 
at a specific galacto-centric radius. There are only two free parameters, the position of the start of 
the flare and the scale length of the flare 

It should be noted that the aim is not to try and fit the model to the data but to demonstrate 
that a simple model does reproduce the counts. Effects such as errors in the extinction calculation, 
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residual effects or the warp and the expected none axi-sysmetry of the outer disc would be require 
a far more complex model and a far larger coverage of sky, in particular closer to the plane. 

In Figures 2 and 3 the effects of various flares starting at a galacto-centric radius of 16 kpc 
are shown. For Fig. 2 the amplitude A has been normalised to the measured counts between m g 
of 16 and 19, which is well well before the flare starts. In Fig. 3, however, the density has been 
normalized to the measured p{r)dr. 

Fig. 4 takes the densities and determines by how much the disc would need to be flared at each 
galacto-centric radius to give the determined density. At the radius of the Sun the relative scale 
height is 1 and it remains constant until a R of 16 kpc, at which point the scale height increases 
rapidly. It is true that the flare is not identical in each case, however they are very similar with a 
typical flare scale length of 4.5 + 1.5 kpc (Fig 4). There are, however, a number of effects which 
have not been taken into account. 

- The warp has been ignored here, but would mean that the average plane position would be 
above b = in the in the / = 150° region and below b = at 220°. In the region 10-18 degrees 
from the plane the counts the flare increases rapidly, and so a small error in the position of the 
plane would significantly reduce the counts at / = 150° b = 15 . 

- The simple model used here assumes a circularly symmetric outer Galaxy. As here we are 
dealing with the thick disc at 6 to 7 scale lengths from the Galactic centre, then if there were a 
10 % change in scale length with Galactic longitude it would lead to a 50% change in counts. 

- The model used here assumes that the parameters of the flare do not depend on Galatic longi- 
tude. 

- The extinction model used here is very simplistic. 

Finally it should be noted that the counts for most distant points are very low leading to a large 
scatter. When far more regions become available it will be possible to quantify these effects and so 
refine the results. 

These results do show that the variation of the number of stars with galacto-centric distance is 
consistent with an exponential law out to over 20 kpc. There is also no evidence for either a major 
drop in counts just before the flare starts or a sudden jump in the counts. This means that the counts 
in the disc can be described by a single function to beyond 20 kpc. 

These plots show that the simple flare model not only predicts that at 20 degrees from the plane 
bf that there should be two peaks in the star counts, but can reproduce the numbers of counts and 
form. 

The alternative scenario is a ring of stars centred at R- 23 kpc, which is detected for R > 16-17 
kpc (Conn et al. 2007). As this structure would only exists in the outer Galaxy, its parameters are 
not constrained to fit observations in the solar neighbourhood or elsewhere in the Galaxy. This 
means that there is a wide range of possible parameter space that can be explored and hence by 
adjusting the width and density, a ring can be made to fit second peak. However, this ring must start 
very close to the edge of the disc. If the inner part of the stream is more than 1 or 2 kpc from the 
edge of the disc then there would be a noticeable drop in the stellar density at that point. Conversely 
if it overlapped with the edge of the disc then there would be a jump in the density. Furthermore, 
the number of stars in the stream has to be very similar to those predicted for the disc with no cut 
off, otherwise none of the flared models would come close to reproducing the measured counts. 
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Finally the vertical distribution of the counts again has to mimic the flared disc model, and this is 
over more that 20 degrees in latitude which correspond to 3 to 4 kpc at an R of 15 to 20 kpc. Hence, 
if there is a stream it would have to have many characteristics which are identical to those predicted 
for a flared outer disc. 

6. Discussion 

That the disc of the Galaxy is flared is widely accepted and is particularly well studied in the HI 
(e.g., Kalberla et al. 2007). They also find that the flare becomes relevant at around R=16 kpc and 
the disc continues to at least 23 kpc. There may be also a dependence of the scale height with the 
galacto centric azimuth, as observed in the gas too (Kalberla et al. 2007), and this lopsidedness 
could find an explanation on the non-axisymmetrical distribution of pressures in the outer disc 
(Lopez-Corredoira & Betancort-Rijo 2009) 

By, in effect, moving sources away from the plane the result is that once the flare starts there 
will be a sharp drop in the counts when compared to the predictions of a simple exponential disc 
for any lines of sight near the plane, whereas at higher latitudes there will be considerably more 
stars than predicted. The sudden loss of stars beyond about /?=14 kpc is reported in Robin et al. 
(1992) when looking at I = 179° b = -2.5°, and the results shown here shows that the second peak 
in the differential star counts at b = 20° is entirely consistent with a flare. 

The stars that are in the second peak would belong to the thick disc as they are typically 3 to 
4 kpc above the plane and are at two to three times the galacto-centric radius of the Sun. Their 
metallicity has been measured to be an average of [Fe/H]=-0.96, with an rms scatter of only ~ 0.15 
dex (Ivezic et al. 2008), which is compatible with the metallicity of the outer thick dick (Lopez- 
Corredoira et al. 2007, Fig. 3). 

The flare of the thick disc is important for understanding its origin, as current theories suggest 
that thick discs caused by mergers will always flare. Bournaud et al. (2009) found that if a flare 
were caused by minor mergers then there would be no region of the thick disc with a constant scale 
height. The analysis presented here is not sensitive to changes in the scaleheight of thick disc closer 
to the solar circle as the counts for R < \6kpc along the lines of site used are dominated by the 
thin disc. An analysis of the thick disc flare in the whole range of galactocentric distances would 
be necessary to see whether the mechanism proposed by Bournaud et al. (2009) is applicable here. 

7. Conclusions 

Uncertainties in our knowledge of the disc makes asserting the presence of the Monoceros stream 
principally because a model does not accurately reproduce the star counts dangerous, even when 
using a model as good as the Besancon model. If the disc has a cut off at R of about 14 kpc, as 
used in the Besancon model, then the ring would be necessary to give the star counts detected. If, 
however, there is no cut off and the thick disc instead flares, then the bulk of the stars in the regions 
being discussed here cannot belong to a stream as there are only sufficient stars to be explained by 
the disc alone. It should be noted that this would not preclude there being smaller scale streams 
towards the anti-centre affecting smaller regions of sky. 

We have shown that a flare in the thick disc would leave a very clear signature in the star counts 
in the outer Galaxy, and that this signature is seen. We have shown that both the form and number 
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Fig. 4. The proportional flare in the scale height required to give the counts at each magnitude along 
each line of sight. The error bars from the number of sources detected in each magnitude bin are 
shown for one region. 



of stars detected towards the anti-centre between 1 1 and 3 1 degrees from the plane can be readily 
explained by a flared thick disc without a cut-off. When a larger area of sky becomes available it 
will be possible to refine the outer disc and flare models and so improve the quality of the fit, in 
particular nearer the plane. Hence, the presence of the extra-galactic Monoceros ring/stream is not 
required to explain the star counts. 
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